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Abstract 

The master equation describing non-equilibrium one-dimensional 
problems like diffusion limited reactions can be written as a euclideen 
Schrodinger equation in which the wave function is the probability 
distribution and the Hamiltonian is that of a quantum chain with 
nearest-neighbour interactions. Since many one-dimensional chains 
are integrable, this opens a new field of applications. For many reac- 
tions the Hamiltonian can be written as the sum of generators of vari- 
ous quotients of the Hecke algebra giving hermitian and non-hermitian 
(for irreversible processes) representations. 



1 Introduction 

It is the aim of this lectures to show that mathematical notions like asso- 
ciative algebras and quantum groups which have appeared in the context 
of integrable systems and equilibrium statistical mechanics find a natural 
application in non-equilibrium statistical mechanics also. It is well known, 
for example, in the literature that several integrable quantum chains corre- 
sponding to magnetic systems [1] can be represented as sums of generators 
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of Hecke algebras H n (q) when certain artificial interactions are added in the 
bulk and at the surface (q is the same parameter which appears in quantum 
groups). We will show that the master equation describing the dynamics of 
some chemical processes limited by diffusion [2] gives representations of Hecke 
algebras, where the supplementary interactions appear naturally and the de- 
formation parameter q occurring in quantum groups gets a simple physical 
meaning. Moreover we would like to stress that the quantities of interest in 
non-equilibrium statistical mechanics (concentrations, correlation functions) 
are very different than those occurring in equilibrium problems and that the 
knowledge obtained in integrable systems to get spectra and wave-functions 
of quantum chains Hamiltonians can find new fields of applications. 

This lectures are organized as follows. In Sec. 2 we review the definitions 
of the Hecke algebras and remind the reader about their known hermitian 
realizations (g-real) through the Perk-Schultz [3] quantum chains. 

In Sec. 3 we describe some reaction-diffusion phenomena and show their 
critical properties. For the case of one-dimensional chemistry (molecules 
diffuse and react on one line - not a very realistic experimental situation) we 
give the master equation and the corresponding Hamiltonian which gives the 
time-evolution of the system. In Sec. 4 we consider the simple example of 
the one-species annihilation model. This example is rich enough to illustrate 
how integrable systems, Hecke algebras and quantum groups appear in non- 
equilibrium phenomena. More examples are given in Sec. 5. The examples 
are of two kinds. Some chemical reactions are given by the same quotients of 
the Hecke algebra as those appearing in the one-species annihilation model 
which implies that the energy levels (not the degeneracies) are the same. 
Other chemical reactions give new quotients of the Hecke algebra. As a 
matter of fact, we show that all the quotients of the Hecke algebra have 
realizations in chemical processes. 

At this point the reader might wonder if the occurrence of Hecke alge- 
bras in reaction-diffusion processes is not just a curiosity, the answer is that 
starting with a realization of the Hecke algebra one can construct through 
"Baxterisation" [4] a spectral dependent matrix R(u) satisfying the Yang- 
Baxter relations which is the first step towards the exact integrability of the 
Hamiltonian. 

These lectures are based on Refs. [5 - 7]. For applications of the ideas 
presented here, see Refs. [8 - 10]. For related topics and complementary 
approaches see Refs. [11 -13]. 
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2 The Hecke algebra and its quotients 



The Hecke algebra H n (q) (with n — L — 1) is an associative algebra with 
generators (i — 1, • • • , L — 1) satisfying the relations: 

eje i± iei - = e i ± 1 e i e i ± 1 - e i± i (2.1) 

[e 4 ,e,]=0, \i-j\>2 (2.2) 

e? = (9 + 9 _1 )ei (2.3) 

where q is in general a complex parameter. In our applications q is always 
real. 

One can define [14] a sequence of (P, M) quotient algebras of H n (q) (P 
and M are non-negative integers). Since the definition is quite complicated, 
we will give a few examples only. For the quotient (1,0), one just take = 0. 
This quotient is not relevant for our purposes. The (2,0) quotient is the 
Temperley-Lieb algebra: 

eie i± iei - e» = (2.4) 

The (1,1) quotient reads: 

(eie i+2 )e i+1 (q + q~ l - e*) (q + q' 1 - e i+2 ) = (2.5) 

Different representations of the (2,0) and (1,1) quotients will be given in 
Sees. 4 and 5. We will not give here the definitions of other quotients but 
just mention that each (P, M) quotient has as representative one of the (P, 
M) Perk-Schultz [3] quantum chains where the Hamiltonian if ( p > M ) is given 
by the sum of H n (q) generators: 

H (P,M) = e f M) (2.6) 



where 



+ q -^ E sign(a - (3) E™Ef +1 ) (2.7) 
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with N = P + M. The N x N matrices E a/3 have elements: 



(E^)^ = 5 ai 5 ps (a, (3 = 0, 1, • • • , N - 1) 



and 



(2.8) 



(2.9) 



£o — e i — ' ' " — £p-i — — £ P — • • • — — Sp+M-l — 1 

The quantum chains ff( p ' M ) are U q (su(P/M)) invariant [15] and have the 
important property that they contain all the irreducible representations of 
the (P, M) quotient. In the examples we are going to mention, if a chain 
corresponds to a certain quotient (P, M) then its energy levels (not the de- 
generacies!) are the same as those of the representative of this quotient given 
by eqs. (2.6) and (2.7). For a detailed discussion of the degeneracies and 
the role of finite-dimensional representations of the affine quantum groups 
U g (su(ri)), see Ref. [7]. 

The generators of the (2, 0) and (1,1) quotients can be written in terms 
of the Pauli matrices a? a y , and a z Aj = 1, • • • , L): 



,(2,0) 



i+i 



2 VjVj+l 



1 fa ~ 



2 



(2.10) 



The Hamiltonian (see eqs. (2.6) and (2.10)) is the U q (su(2)) invari- 

ant quantum chain of Pasquier and Saleur [16]. The U q (su(2)) generators: 



1 



/ 3=1 



± 



3=1 
IV 



IV' 

2 Z^A 



commute with iP^ 2 '°) and satisfy the relations 

= [2s 2 ], 



q 2S* _ q -2S* 



1-1 



(2.11) 

(2.12) 
(2.13) 

(2.14) 
(2.15) 
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These Hecke generators of the (1,1) quotient are: 



1 



ef' 1} = "2 [ ° X j°U + <#3+i + + K " (9 + <T 1 ) ] (2-16) 

The Hamiltonian H^ 1 ' 1 ^ commutes [17] with the generators of the super- 
algebra U q (su(l/l)),S z (see eq. (2.11)) and T ± : 



3=1 



S Z ,T± 



{T + ,T-} 



q L - q~ L 



q-q- 



(2.17) 
(2.18) 



That much for the mathematical introduction. 



3 Reaction - diffusion phenomena 

One of the simplest examples of a reaction-diffusion process is the one in 
which one has two types of molecules, say A and B which diffuse in a milieu 
(a gel solution) and react to give an inert product which precipitates. If at 
the time t — the concentrations of the molecules A and B are 

C A (t = 0) = C B (t = 0) = ip A (3.1) 

one observes that for long times, one has: 

C A (t) = C B (t) = A (3.2) 

where d is the space-dimensions of the milieu and A is a constant independent 
on the initial concentration (p A - The independence of the late-times concen- 
trations on the initial conditions is called self-organisaion. Notice the power 
fall-off of the concentrations which implies a critical (massless) behaviour. 
For the experimental setting and results see Ref. [18]. There are many theo- 
retical and experimental questions to ask about such a process, we are going 
to mention a few at the end of this Section. 

In order to model the process just described in d — 1, we will consider a 
one-dimentional chain with L sites. To shorten the notations we will denote 
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by Aq a vacancy and by A\ and A2 the molecules A and B respectively. At 
t = we give a probability to find a distribution of the molecules Ap(f3 = 

0. 1.2) on the L sites (a vacancy A will also be called a molecule). In 
order to describe the time evolution of the system we will assume that there 
are only two-bodies nearest-neighbours "interactions" decribed by the rates 
pa,p which give the probability per unit time for two molecules A a and Ap 
situated on the sites i and i + — 1, - • • ,L — lto change into the molecules 
A 1 and A a . Notice that in general the processes don't have to be left-right 
symmetric (r"'^ 7^ ■ F° r ^ ne process described above one can take only 
a few non-zero rates: 

A k + A -> A + A k (diffusion to the right) (k = 1, 2) (3.3) 

A + A k -> A fc + A (diffusion to the left) (k = 1, 2) (3.4) 

Ai + A 2 -> A + Ad (reaction) (3.5) 

A + A -> A + A (reaction) (3.6) 

with only two independent rates D and Y. 

r fc ; ° = r° fc ; fc = d (k = 1, 2) , r^;^ = rj|;J = r (3.7) 

The time evolution of the system is given by a master equation which 
gives the probability to find a certain configuration of the molecules Ap at a 
time t on the L sites for a given probability distribution at t — 0. Knowing 
the probability distribution at t, one can compute various average quantities 
like the concentrations for L sites and then take the large L limit in order to 
get the result given by eq. (3.2). At this point we will consider the problem 
in its full generality. 

Let us consider an open chain with L sites, where at each site i = 

1, 2, • • • , L we attach a variable (5 (representing the molecule Ap) taking N 
integer values ((3 = 0,1, ■ ■ ■ , N — 1). By convention we attach the value (3 = 
to a vacancy (inert state). The master equation describing the time evolution 
of the probability destribution 



is: 

dP{{(3},t) 



dt 



P({f3};t) = P(f3 1 ,{3 2 ,---,p L ;t) (3.8) 
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N-l 

+ Yl V^i.y^ftrlj+lk. \P j+1 +m] N ,---,P$Q) 
l,m=0 

where the w^ m are related to the transition rates F"'P through the relations: 

T«f = w^ a ^(a,(3) (3.10) 

In eq. (3.9) the prime in the second sum indicates the exclusion of the pair 
I — m — and the symbol [x+2/]tv means the addition (x+y) modulo N. The 
rate u>oo(ct, /3) gives the probability per unit of time that the configuration of 
two neighbouring sites (cx,f3) is changed. Using eq. (3.10), the conservation 
of probabilities implies: 

w , (a,P) = J2w r , s (a - r,0 - s) (3.11) 

r,s 

where r = s = is again excluded. The notations used in eq. (3.9) are 
convenient to read symmetries. The rates are always non-negative and sym- 
metries imply that either some rates are equal or that they vanish. For 
example to have Zn symmetry implies that all rates wi jm (a, (3) vanish except 
those satisfying the relation: 

[l + m] N = (3.12) 
Parity conservation (left - right symmetry) implies: 

wi, m (a, f3) = w m j(P, a) (3.13) 
The master equation (3.9) can be interpreted as a Schrodinger equation 

d 



dt 



P>= -H\P> (3.14) 



in Euclidean time if we identify the probability distribution P({(3}; t) as the 
wave function. The Hamiltonian in eq. (3.14) 

H=J2H J (3.15) 

3=1 
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acts in a Hilbert space of dimension N L while H 3 acts in the subspace <g> 
V^ +lS) and is given by 

H 3 = U 3 -T 3 (3.16) 

where 

N-l N-l 

T j = E' E Wi,m(a,/3)^ ,[a+,]jv ®E^ +m ^ 

l,m=0 a,/3=0 
N-l 

Uj = E w , (a,P)E aa ®EW (3.17) 

a,/3=0 

The matrices E a @ have been defined in eq. (2.8). We now show how to 
compute physical quantities from the knowledge of | P >. We define a basis 

| {(3} >=\ p u fo, ■ ■ -,p L > , < {(3'} | {(3} >= 5 mm (3.18) 

and is this basis the ket-vector | P > is: 

| P >=$>({/?},*) !{/?}> (3.19) 

{/3} 

to the initial probability distribution P({{3},t = 0) corresponds the ket- 
vector | P >. 

We define a left-vacuum bra-state: 

< |= E < {0} I (3-20) 

{/3} 

then the conservation of probabilities (3.11) implies 

<0\H = (3.21) 
If X({/3}) is an observable, then its average at a time t is 

< X > t = X({P})P({Ph t)=<0\ Xe~ m I P > (3.22) 

{/3} 

Notice that the expression (3.22) is very different of the vacuum expecta- 
tion values encountered in equilibrium statistical mechanics. Here the initial 
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probability distribution | Pq > appears in an essential way and can introduce 
new lenghts scales into the problem. 

If the energy levels E x and the eigen-functions | ip\ > are known 

H | ip\ >= E x \ip\> (3.23) 

with 

I P >=E«a |^a> (3.24) 

A 

< X > t can be rewritten as: 

<X > t =J2 a\e~ Ext < | X | ^ A > (3.25) 

A 

The system has a critical behaviour if in the thermodynamic limit (L — > oo) 
all the E\ vanish (keep in mind that because of eq. (3.21), the ground state 
has always energy zero). There is a profound difference between a massless 
regime often (but not always!) seen in equilibrium statistical mechanics when 

lim 

LE\ = constant (3.26) 

L — > oo 

and one has conformal invariance and the present case when 

r hm L 2 E X = constant (3.27) 
L — > oo 

This can be easily understood since the diffusion equation looks like a non- 
relativistic Schrodinger equation and eq. (3.27) reflects a non-relativistic 
(quadratic) dispersion relation. Quadratic dispersion relations are met also 
in equilibrium problems like in the Pokrovsky-Talapov [19] phase transition. 

Before closing this Section, let us observe that it is often very useful 
to perform similarity transformations in the master equation (3.9). This 
implies considering instead of P({(3};t), another function ijj({(3},t) related 
to the first through the relation: 

P({(3};t) = ^m,t) (3.28) 

where $ acts in the {f3} space but is time independent. For our purposes 
here, we take $ diagonal and of a factorized form: 

^=f[h (k \p k ) (3.29) 
k=i 
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The new master equaiton reads: 



dt ~ ^= 1V " ,UM '' '■' 



Zi=a-woM,Pj+iM{P},t) + 



where 



Accordingly, after the similarity transformation, the new Hamiltonian (3.15) 
- (3.16), has the wi m (a, (3) in the expression of the 7}'s replaced by the W^s 
given by eq. (3.31). 



4 The annihilation-diffusion process 

We now present a simple example in which a single spades of molecules A\ 
annihilate: 

A 1 + A 1 -> A + A (4.1) 
with a rate u>i 5 i(0, 0) = a, diffuse to the right 

A 1 + A -> A + Ax (4.2) 

with a rate u>n(0, 1) = Dr, and diffuse to the left 

A + A x -> A + A (4.3) 

with a rate w n (l,0) = D^. No other processes are considered. We will 
denote 

D = ^D L D R , 9 = ^ (4-4) 
and make a similarity transformation (see eq. (3.29)): 

fr (fc) (l) _ 9^ ^ 

fc(*)(0) A 1 J 
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where A is a free parameter. We also choose the time units such that D = 1 
such that the whole process depends on two physical parameters a and q. 
The Hamiltonian describing this process is: 

H = H + H 1 (4.6) 

where 

L-\ 

2 



where 



- ^ - - + 2A' -A - 2] (4.7) 

A = ^-f, A' = l-| (A'<A) (4.9) 

Notice that H is hermitian and that the non-hermitian H 1 has a coupling 
constant proportional to A~ 2 and that the parameter A is given by the simi- 
larity transformation (4.5). Since the spectrum of the Hamiltonian H can't 
depend on a similarity transformation this implies that the spectrum of H 
coincides with the spectrum of H . The spectrum of H is known since this 
is the XX Z chain in an external Z field [20] where the Bethe-Ansatz applies. 
Leaving asside the surface fields, H can be written as 



K = ~\ E H + + A.] (7? +1 + Ka] + (4.10) 

z j=i 

In the thermodynamic limit (L — > oo) the spectrum of i^o is massive if 
h — 1 — A' > 1 — A. This implies that for assymmetric diffusion (q ^ 1) in 
an open chain the spectrum is always massive. If h — 1 — A (which implies 
q — 1) the spectrum is massless and corresponds to a Pokrovsky-Talapov 
phase transition with a quadratic dispersion relation. 

It is now time to make the connexion between the annihilation-diffusion 
process and the Hecke algebra. Let us first consider pure diffusion (a = 0). 
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It is easy to check that in this case (see eq. (2.10)) that 

# = tfo = Ee?' 0) (4.11) 
3=1 

This result is very interesting since it gives a simple physical interpretation 
to the parameter q of the quantum group U q (su(2)). It is (see eq. (4.4)) 
the square root of the ratio of the left and right diffusion constants. The 
spectrum is of course massive if q ^ 1. This is true for the open chain. If 
one considers the assymmetric diffusion process on a ring however (periodic 
boundary conditions), the similarity transformation (4.5) induces a boundary 
term with a strength of the order of the volume which makes the system 
massless [21]. 

Let us now take 

a = q + q~ l (4.12) 

which corresponds to the physical situation if two molecules A 1 are on neigh- 
bouring sites, they annihilate with probability one. In this case 

Ho = E ef» (4.13) 

3=1 

moreover 

H = E fj (4.14) 

3=1 

where the non-hermitian matrices fj satisfy the relations (2.1 - 3) and (2.5) 
giving a representation of the (1,1) quotient of the Hecke algebra. 

The example considered in this Section illustrates the general structure 
which is going to be described in the next section. Pure diffusion leads to 
quotients of the type (N, 0) corresponding to U q (su(N,0)) symmetric quan- 
tum chains, non-reversible processes give (M, P) (P ^ 0) quotients. 

5 Various representations of the Hecke alge- 
bra given by reaction-diffusion processes 

(a) Diffusion processes. We consider (JV— 1) types of molecules (A±, ■ ■ ■ A^-i) 
which diffuse to the right: 

A b + A ^ Ao + A b (&=!,-• -,N- 1) (5.1) 
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with equal rates Tr and to the left: 

A + A b ^A h + A (b = l,---,N-l), (5.2) 

with equal rates T L . Choosing the time scale in such a way that the diffusion 
constant 

D = y/TiTn = 1 and introducing q = ^JTr/Tl, up to an equivalence trans- 
formation, we obtain Hj = ej corresponding to the quotient (2,0). 

(b) Diffusion and interchange processes. To the diffusion processes (5.1) 
and (5.2) we add the interchange to the right processes: 

A b + A C ^A C + A b (b^c,b>c), (5.3) 

with rates T R and with rates T L if b < c. In this case the quotient is (N, 0) 
and the obtained representation of H n (q) is hermitian (the same applies to 
the diffusion processes of 

(c) Diffusion, interchange and coagulation processes. We now add to the 
processes (5.1 - 3) only one reaction for the molecule A 1 (the molecules are 
ordered!): 

A 1 + A 1 ^A + A 1 , (5.4) 

with a rate T R and the process where A and Ai interchange their position 
in the final state with rate I\. In this case the quotient is (N — 1,1) and the 
representation of H n is non-hermitian. As a rule, non-hermitian representa- 
tions occur for quotients (P, M) with M ^ i. e. when the representative 
Hamiltonians if ( p ' M ) are supersymmetric. 

One obtains new quotients of H n generalizing the reaction (5.4) by taking 

A r + A r -> A s + A r (ratelfl), (5.5) 

A r j + A r j -> A 8 , + A r , (rater^), (5.6) 

where s — r ± 1, s' — r' ± 1, and r ^ r', s'; s ^ r', s'. In this way one obtains 
the quotient (N — 2,2). Taking more reactions like in eqs. (5.5), (5.6) one 
can get arbitrary quotients (P,M). 

(d) Diffusion-polymerisation (modulo N) processes. In this case besides 
the diffusion (eqs. 5.1), (5.2), one considers the processes 

A r + A s — > Aq + A[ r+S ] N 
(rater R ,[r + s] N ^0), (5.7) 
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(rate r L , [r + s] N ^ 0), 



(5.8) 



A + A s ^A + A 
(rate + F L , [r + s] N = 0) 

(r,s = l,2,...,JV-l). (5.9) 

For all these processes one obtains the quotients (1,1). 

(e) Diffusion, interchange and annihilation processes. Because of inter- 
change processes, the labels of the molecules are ordered; we can now add 
the following two processes: 

A m + A m -> A m+1 + A m+1 (rateAi), (5.10) 

A m + A m -> A m _i + A m _i (rateA_i), (5.11) 

with Ai + A_! = + T^. In this way we obtain again the quotient (2, 1). 

The reader might wonder how the above H n (q) quotients were obtained. 
From the work done in ref. [5] we have developed a certain intuition for 
chemical processes up to A^ = 3. We have then used the computer to check 
that no other solutions exist for N = 2 and 3 and generalized the solutions 
obtained for A^ up to 3 to all A^ and checked them. There are certainly other 
solutions which one can discover studying sytems with A^ = 4 and certainly 
examples of the quotients (2,2) different from the one given in eqs. (5.5), 
(5.6) will show up. We hope that the present effort to connect representatives 
of H n (q) and chemical processes will help to compute various correlation 
functions for the latter. One immediate consequence of the above mentioned 
connection is that all processes with g ^ 1 are massive (the concentrations 
for example have an exponential fall-off in time) and all those with q — 1 are 
massless (algebraic fall-off). 
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